Nonlinear Relativistic and Quantum Equations with a Common Type of Solution 
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Generalizations of the three main equations of quantum physics, namely, the Schrodinger, Klein- 
Gordon, and Dirac equations, are proposed. Nonlinear terms, characterized by exponents depending 
on an index q, are considered in such a way that the standard linear equations are recovered in the 
limit q — >■ 1. Interestingly, these equations present a common, sohton-like, travelling solution, 
which is written in terms of the g-exponential function that naturally emerges within nonextensive 
statistical mechanics. In all cases, the well-known Einstein energy-momentum relation is preserved 
for arbitrary values of q. 
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The applicability of linear equations in physics is usu- 
ally restricted to idealized systems, being valid for me- 
dia characterized by specific conditions, like homogene- 
ity, isotropy, and translational invariance, with particles 
interacting through short-range forces and with a dynam- 
ical behavior characterized by short-time memory. How- 
ever, many real systems - specially the ones within the 
realm of complex systems ~ do not fulfill these require- 
ments, and usually exhibit complicated collective behav- 
ior associated to NL phenomena [ij. Accordingly, the 
study of nonlinear equations (NLEs) has opened a new 
area of physics, attracting a lot of interest due to the pos- 
sibility of describing many real situations. Since finding 
analytical solutions of NLEs can be a hard task^partic- 
ularly in the case of NL differential equations [2|, very 
frequently one has to make use of numerical procedures 
and so, a considerable advance has been attained lately. 

Among the most studied NL differential equations we 
have the sine-Gordon the NL Schrodinger Q and 
Fokkcr-Planck ones Q. In these cases, the NL contribu- 
tions lead to important novel insights, relevant for mod- 
eling several physical new aspects. Two special types 
of solutions have called a lot of interest, particularly in 
the sine-Gordon and NL Schrodinger equations; these are 
known as solitons (or solitary waves) and breathers [ll-Q ■ 
In both cases one has compact travelling solutions, the 
first one being characterized by a spatial part that does 
not deform throughout the evolution, whereas the later 
presents an oscillating width as propagation occurs. An 
important property of these two solutions concerns their 
square integrability, allowing for an appropriate normal- 
ization. Due to the modulation of the wave function, 
these solutions are considered to be relevant in diverse 
areas of physics, including nonlinear optics, supercon- 
ductivity, plasma physics, and deep water waves 

In the most common NL Schrodinger equation, one 



introduces a new cubic term in the wave function (see 
for instance Q ) which, for some particular type of solu- 
tion, is responsible for the modulation of the wave func- 
tion. Such addition of extra NL terms has been much 
used in the literature for constructing NLEs. A differ- 
ent approach consists in transforming one or more linear 
terms into NL ones, as usually happens in NL Fokker- 
Planck equations [5[. In this case, whereas the linear 
Fokker-Planck equation is associated to normal diffusion 
and to the Boltzmann-Gibbs entropy, its NL counterparts 
are usually related to anomalous-diffusion phenomena @ 
and to generalized entropies , fike the nonadditive one 
which yields nonextensive statistical mechanics 

This theory emerged from a generalization of the 
Boltzmann-Gibbs entropy, through the introduction of 
an index q {q 1 recovers the standard case). Since 
then, a considerable progress was achieved, leading to 
generalized functions, distributions, important equations 
of physics, and even new forms of the Central Limit 
Theorem 11|. In particular, the g-Gaussian distribu- 



tion, which generalizes the standard Gaussian, appears 
naturally by extremizing the entropy Q, or from the 
solution of the corresponding nonlinear Fokker-Planck 
equation [l2l |. and has been very useful for experiments 
in many real systems [l3|- Among many others, we 
may mention: (i) The velocities of cold atoms in dissi- 
pative optical lattices ; (ii) The velocities of particles 
in quasi- two dimensional dusty plasma (iii) Single 
ions in radio freq uency traps interacting with a classi- 
cal buffer gas [15| ; (iv) The relaxation curves of RKKY 
spin glasses, like CuMn and AuFe [l^; (v) Transverse 
momenta distributions at LHC experiments [l7j . 

Herein we introduce NL generalizations of the three 
main equations of quantum physics, namely, the 
Schrodinger, Klein-Gordon, and Dirac equations. The 
present proposals consist in extending linear terms into 
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NL ones, as done for the NL Fokkcr-Planck equations. 
An interesting aspect about these generahzations is that 
they may be formulated easily in arbitrary dimensions, 
whose exact solutions are here presented. These solutions 
are all expressed in terms of the g-exponential function 
,(u)th 

Urn. 

as the principal value of 



exp^(M)that naturally emerges in nonextensive statis- 
tics ISlllOl- For a pure imaginary iu, one defines expg(Mt 



exp^(m) = [1 + (1 — q)iu]~ ; expj(m) = exp(m). (1) 
The above function satisfies [l8| . 

exp^(±iM) = coSqiu) ± 2sing(u) , (2) 
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coSg(u) = Pq{u) COS < j-arctan[(q — \)u\ 



(3) 



sinq(7i) ~ Pq{u) siu ^arctan[(g — \)u] ^ , (4) 



(5) 



exp^(m) expg(-m) = cos^(u) + sin^(M) = . (6) 

Notice that exp^ + U2)] 7^ exp^(mi) exp^(m2) for 
<7 7^ 1 As a consequence of Eqs. the q- 

exponential of a pure imaginary presents an oscillatory 
behavior with a varying amplitude Pq{u) that decreases 
(increases) for 1 < g < 3 (g < 1). By integrating Eq. (|6]) 
from —00 to +00, we verify the physically important 
property of square integrability for 1 < q < 3, whereas 
such integral diverges in both limits g — 1 and g -> 3 
and for g < 1 [l^ . More precisely, exp^ (m) is modulated 
by the g-Gaussian in Eq. ([5]). This is the type of solution 
that we focus on throughout the present paper. 

Let us first consider the simple one-dimensional linear 
wave equation 



(7) 



for which any function of the type — ujt), twice 

differentiable, is a solution. In particular, one may have 
a q-plane wave [l9l |. 

$(a;,t) = $0 expg[i(fcx-Lj<)] ; [$o = $(0,0)], (8) 

as a solution of the equation above, provided that co ~ 
ck. Let us remind that the above g-plane wave cannot be 
decomposed, for g 7^ 1, into spatial and ternporal factors, 
as it happens in many solutions of NLEs [2|. Moreover, 
the amplitude of the g-plane wave for 1 < g < 3 de- 
creases when its argument {kx — ujt) increases. Due to 
this property, this new type of solution may well be ap- 
plicable to various nonlinear physical phenomena. Let us 
stress that, in the particular situation where x = ct, one 



has ^{x,t) = $0 (Vt), consequently the g-planc wave be- 
haves like a soliton propagating with a velocity c = uj/k. 
This enables the approach of nonlinear excitations which 
do not deform in time and should be relevant, e.g., in 
nonlinear optics and plasma physics. 

Now, since Eq. ([7]) is a linear equation, the following 
linear combination. 



$(a;, t)=Y^ $oj cxpg^ [i{kx - ujt)] 



(9) 



is also a solution with the same dispersion relation. A 
given physical system may be characterized by a single 
value of g. If 1 < g < 3, we have ^{±00, t) = (Vt); by 
specifying the rate of decay of the g-plane wave ampli- 
tude (which should be a characteristic of a given physical 
system), i.e., its modulation, Eq. ([5|), we may determine 
the appropriate value of g. 

Another important aspect of the g-plane wave concerns 
its immediate extension to d dimensions. 



$(a?, t) = $0 exPg i{k ■ X — ut) 



(10) 



which, due to the well-known properties of the scalar 
product k ■ X = X]n=i^n-^"i exhibits invariancc under 
rotations. If we take into account that dcxp^{z)/dz = 
[expg(z)]'^ and d'^ expy{z)/dz'^ = q[cxp^{z)]'^''^^ wc ob- 
tain, for the d-dimensional Laplacian, 



VH{x,t) = -q (^Zj^nj $o{cxpJz(fc ■ 



x-ujt)]} 



2q-l 



(11) 

^From the results above one sees that the d-dimensional 
g-plane wave of Eq. (fTOl) satisfies the linear wave equa- 
tion. 



V2$(x,t) 



1 d^<i>{x,t) 



ck 



(12) 



with k = Y X]n=i 1 showing that the resulting disper- 
sion relation is g-invariant. 

We will now focus on NLEs for which the above d- 
dimensional g-plane wave is an exact solution. It should 
be emphasized at this point that many NLEs in the liter- 
ature are usually formulated in one dimension and that 
their extension to d dimensions may not be always an 
easy task. Accordingly, let us introduce the following d- 
dimcnsional NL generalization of the Schrodinger equa- 
tion for a particle of mass m. 



d 
'Ft 



$(x,t) 



1 



2 - g 2m 



<^{x,t) 



l-2-q 



(13) 



We notice that the scaling of the wave function by $0 
guarantees the correct physical dimensionalities for all 
terms. This scaling becomes irrelevant only for linear 
equations [e.g., in the particular case g = 1 of Eq. p3p ]. 
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Consistently, the energy and momentum operators are 
generalized a.s E = iTiDt and p„ — —ifiDx^ respectively, 
where Duf{u) = [f[u)Y^''df{u)/du. These operators, 



when acting on the g-exponcntial exp^ 



i{k ■ X — uit) 



yield the energy E — fiuj and momentum p = fik. Now, 
if ones considers the g-planc wave solution of Eq. (fTO|) 
by using k — > p/fi and lj — !> E/h^ one verifies that 
this new form is a solution of the equation above, with 
E = /2ra, for all values of q. The NL Schrodingcr 
equation of Eq. ([T^]) shows the same structure of the NL 
Fokker-Planck equation of Rcfs. [l2| in the absence of an 



external potential, which appears in nonextensive sta- 
tistical mechanics [1, [^. Essentially, it represents the 
analogue of the porous-medium equation l20| , very com- 
mon in anomalous-diffusion phenomena [6[, but with an 
imaginary time. One well-known solution of the porous- 
medium equation is due to Barenblatt 2l| . This solution 



presents a structure very similar to a q-Gaussian, written 
in terms of a power law like the one in Eq. ([8|). 

However, Eq. p3p differs from previous formula- 
tions Q, where one adds a new NL term (in most cases, 
a cubic nonlinearity in the wave function) to the two 
existing linear terms. This extra NL term is responsi- 
ble for the modulation of the wave function, which has 
made such equations relevant for many physical phenom- 
ena since some decades. Particularly, the breather-like 
solutions are characterized by a time dependence of the 
linear type, i.e., e*"*, such that their nonlinearities are 
essentially manifested through their spatial dependences. 
In its discrete formulation, this equation leads to discrete 
breathers, much studied in the literature due to their po- 
tential applications in lattice dynamics of solids, coupled 
arrays of Josephson junctions, and localization of waves 
in random lattices, among others (22I. l23l] . The main dif- 
ferences between Eq. (fT3|) and the existing forms of NL 
Schrodinger equations in the literature are: (i) Instead of 
adding an extra term in which the nonlinearity is intro- 
duced, we modify the spatial second-derivative term; (ii) 
The equation, together with the proposed solution, are 
easily extended from one to d dimensions; (iii) The cor- 
responding solution of Eq. ()13|) manifests nonlinearity in 
both space and time, through a modulation in these two 
variables, which keeps the norm finite for all {x,t); (iv) 
The well-known energy spectrum E = p'^ /2m is preserved 
for all q. Therefore Eq. ([T3| . together with the simple 
solution in Eq. ([T0|) . emerge as possible descriptions of 
physical phenomena, like the propagation of nonlinear 
pulses in optical fibers for carrying information ( "bits" ) . 

Another category of NL equations studied in the liter- 
ature concerns the class of NL Klein-Gordon equations 
(see for instance (il. Isl. [23l - [25t ) . As in the previous case, in 
most NL Klein-Gordon equations the second-derivative 
terms are left unchanged and the nonlinearity is intro- 
duced by means of extra terms containing powers of the 
Klein-Gordon field e.g., cubic or quartic 



terms. If one extends this nonlinearity to a general func- 
tional of the field, then the sine- Gordon equation 
which contains a term expressed as a sine of the field, 
may be also included in this category [l|, Several 
types of solutions have been proposed in the literature 
for such NL Klein- Gordon equations, particularly those 
of the Barenblatt kind 21 1 and breather- like ones 22 , 2^ . 
In other cases the solutions are characterized by an ampli- 
tude that may decay in both space and time, multiplied 
by a periodically oscillating part in the form Q^i^^+^t) 
{A and B constants) UM- The solution presented 
in [24, 25 1 is qualitatively similar to the one-dimensional 
q-plane wave. However, it may not be trivially extended 
to d-dimensions. 

Herein we propose a new NL Klein-Gordon equation 
in d dimensions, namely 



V2$(f,t) 



1 d^^{x,t) 
c2 dt^ 



$(f,t) 



2(9-1) 



(14) 



One may verify easily that the same g-plane wave used 
for the NL Schrodinger equation is a solution of Eq. (|14[) , 

preserving for all q the Einstein relation 



E' 



2 2 I 2 4 

p c + m c 



(15) 



It should be mentioned that, to our knowledge, the non- 
linear term in Eq. (jl4p is new and different from those of 
previous formulations, even in the particular cases q = 2 
and q = 5/2, where one gets cubic and quartic terms 
in the field, respectively. In previous works, the nonlin- 
ear term is constructed by multiplying the wave func- 
tion by a power of its modulus, leading to different types 
of solutions and different energy spectra. Another im- 
portant aspect of Eq. (|14p concerns its Lorcntz invari- 
ancc: since this property is directly related to the first 
two terms (containing derivatives) 27[, which have not 
been changed herein, then Eq. (|14p remains invariant un- 
der Lorentz transformation. Some previous NL Klein- 
Gordon equations are not invariant under usual Lorentz 
transformation 



I25| and require its generalization [25j. 
Along the same line we introduce a NL form for the 
Dirac equation. The present proposal represents, to the 
best of our knowledge, a new one (see Ref. [26| for a 
previous, different formulation). In this case, we will re- 
strict ourselves to d ^ 3 spatial dimensions; let us then 
introduce the following generalized Dirac equation, 

in ^^^^ +inc(d-V)<^(x,t) = l3mc^A'^i\x,t) $(f,i) , 
ot 

. 

where ax,ay, ocz (written in terms of the Pauli spin ma- 
trices) and /3 (written in terms of the 2x2 identity ma- 
trix /) are the standard 4x4 matrices [27[. The new, 
g-dependent, term is given by the 4x4 diagonal matrix 
A|j^(x,i) = bij^j(x,t^laj\'^~^, where {oj} are complex 
constants i^J^^^ (x^t) = Sij). The solution of Eq. ()16p we 
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focus on is the following four-component column matrix 



$(f,0 = 



exp„ 



-{p-x-Et) 



of this theory are still elusive. 

We thank Prof. E. M. F. Curado for fruitful conver- 
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(17) 
la, we 



Substituting this four-component vector into Eq. 
get, for the coefficients {a^}, precisely the same set of four 
algebraic equations corresponding to the linear case (see 
page 803, Eq. (15.45b) of (l^). These equations have, 
for all g, a nontrivial solution only if the Einstein energy- 
momentum relation (Eq. ()15|)) is satisfied. The proposal 
in Eq. (|T6|) differs from that in [51] in the same sense of 
the herein proposed NL Schrodinger and Klein-Gordon 
equations, i.e., the nonlinearity is introduced by general- 
izing an existing term, rather than adding an extra NL 
term. As a consequence, our solutions are given by the 
above-mentioned ^-exponential, thus preserving relation 
(fT5|) , in contrast to those in Ref . [2^ , which are typically 
written in terms of the standard exponential. The pos- 
itive and negative parts of the present energy spectrum 
are naturally expected to be respectively associated to a 
particle and its corresponding antiparticle. Contrary to 
the standard plane waves (case q — 1), one has g-plane 
wave solutions, which are square- integrable, i.e., with a 
finite norm for 1 < q < representing localized particles. 

To conclude, we have shown that the g-plane waves, 
which consist of a generalization - within the framework 
of nonextensive statistical mechanics ~ of the standard 
plane waves, represent a very general type of solutions 
for some physical nonlinear equations. Besides the d- 
dimensional linear wave equation, it appears also as a 
solution of important nonlinear generalizations of the 
main equations of quantum physics, namely, nonlinear 
d-dimensional Schrodinger and Klein-Gordon equations, 
and a three-dimensional Dirac equation (its generaliza- 
tion to d-dimensions should in principle be straightfor- 
ward; see [2^ for the linear case). In contrast to other 
generalizations known in the literature, where nonlinear 
terms are added to the usual linear ones ~ and whose 
extensions from one- to d-dimensional space are in many 
cases not particularly straightforward - the present pro- 
posals may be easily formulated in d dimensions. They 
consist in fact in modifying existing terms by introduc- 
ing powers characterized by an index q. The standard 
linear equations are recovered in the limit g — > 1. Fur- 
thermore, the equations and solutions presented here pre- 
serve, in all cases, the usual energy-momentum relations. 
Due to its simplicity and properties described herein, this 
type of solution represents a good candidate for describ- 
ing several nonlinear physical phenomena characterized 
by oscillatory motion with modulation in both space and 
time, like those appearing in superconductivity, plasma 
physics, nonlinear optics, and lattice dynamics of solids. 
Naturally, at the present early stage, specific applications 
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